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SECTION 1 


INTRODUCTION 

This User's Manual describes the input and output variables as well as the Job 
Control Language necessary to utilize the IMP-H apogee motor firing program, 
IMPMOT. The IMPMOT Program can be executed as either a stand-alone pro- 
gram or as a member of the Flight Dynamics System (FDS). This program is 
used to determine the time and attitude at which to fire the IMP-H apogee 
boost motor. 

The IMPMOT Program is written in FORTRAN IV for use on the IBM 360 
series computer. It was created by modifying the existing RAEMOT Program 
(See Reference 1) in the following ways: 

1. The branches necessary to bypass the circular orbit targeting logic 
were created. The alternate initial velocity guidance which was 
implemented is described in Appendix A. 

2. For each apogee motor ignition time, the logic necessary to perform 
a parametric scan upon firing attitude was developed. 

3. The interfaces with the FDS were established. 

4. An 8th order Runga-Kutta integration subroutine was implemented in 
order to propagate the satellite state vector through the apogee motor 
burn. 

5. The logic to accommodate a multisegment apogee motor thrust pro- 
file was developed. 
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6. The logic necessary to create and maintain the sequential excursion 
file was developed. This file contains the satellite state vectors after 
the apogee motor burn for each combination of ignition time and firing 
attitude examined. 

The IMPMOT Program has two guidance schemes available at the users re- 
quest. 

1. Circular orbit guidance 

2. Initial velocity guidance 

1. 1 Circular Orbit Guidance Technique 

When using the circular orbit guidance technique, the program will achieve a 
circular orbit by trading off orbit inclination and nodal position for eccentric- 
ity. Due to the fixed amount of delta-V supplied by the solid propellant apogee 
motor, one of three firing modes will exist. 

1. Excess energy— in which the delta-V supplied by the apogee motor 
exceeds the delta-V required to achieve a circular orbit. 

2. Sufficient energy— in which the delta-V supplied by the apogee motor 
equals the delta-V required to achieve a circular orbit. 

3. Deficient energy— in which the apogee motor does not supply enough 
delta-V to achieve a circular orbit. 

For the excess energy mode, two firing attitudes are possible. The IMPMOT 
Program computes the in-plane delta-V necessary to establish a circular 
orbit and dissipates the excess apogee motor delta-V out-of-plane (hence, the 



change in orbit inclination and node). This out-of-plane component can be 
directed either above or below the transfer trajectory plane (which explains the 
two possible firing attitudes). The program will examine each of the two 
firing attitudes. 

For the sufficient energy mode, only one firing attitude is possible. The pro- 
gram will apply all of the apogee motor delta-V such that circular orbit veloc- 
ity can be obtained. 

The deficient energy mode indicates that a circular orbit can no longer be 
obtained. 

Greater detail of the circular orbit guidance technique can be found in the 
RAEMOT documentation in Reference 1. 

1. 2 Initial Velocity Guidance Technique 

When using the initial velocity guidance scheme, the program maximizes the 
magnitude of the in-plane velocity vector by firing the apogee motor along the 
velocity vector at the ignition time. Appendix A of this document describes 
this technique in detail. 

1. 3 Program Capabilities 

Regardless of the guidance scheme employed, the user has the capability to 
exercise the following options. 

1. A scan upon apogee motor ignition time, tj., can be performed in the 
range; 
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where: 


is the time at which the satellite reaches apogee of the transfer 
trajectory 

Atg is input when using the initial velocity guidance scheme; or is 
determined based upon the delta-V supplied by the apogee motor 
when using the circular orbit guidance scheme.* 

2. For each ignition time defined above, the user may elect to fire the 
apogee motor in a fixed direction. This can be accomplished by in- 
putting the desired firing attitude. Whenever the firing attitude is 
input, the selected guidance scheme is overridden. 

3. The user may elect to examine only one ignition time. This can be 
accomplished by inputting the desired firing time. In this case the 
ignition time scan (defined in item 1 above) is not performed. 

4. For each combination of ignition time and attitude (either computed 
or input), the effects of attitude uncertainties upon the orbit resulting 
from the apogee motor bum will be examined by superimposing the 
attitude uncertainties upon the firing attitude. The following attitude 
combinations are examined for each ignition time 


*When using the circular orbit guidance scheme, the ignition time scan is 
terminated whenever the delta-V supplied by the apogee motor is insufficient to 
establish a circular orbit. 
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where 






0!p. + U , 5p, 
U ot ’ D 


“d 

»D .«D+U, 


“d .«o-U. 


is the spin axis right ascension** 
is the spin axis declination 
is the uncertainty in spin axis right ascension 
Uj is the uncertainty in spin axis declination 
5. In addition to examining the effects of attitude uncertainties, the 
user can perform a parametric scan upon firing attitude about the 
input or computed firing attitude at each ignition time. The range of 
the attitude excursion is defined by: 

«D - 

- ma^ <5<6p +ma, 

where: 

is the increment for the excursion upon right ascension 

a, is the increment for the excursion upon declination 
6 


**The IMPMOT Program assumes that the thrust vector is coincident with the 
negative spin axis. The apogee motor firing attitude is specified by the spin 
axis right ascension and declination. 
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n and m are arbitrary input integers. 

In this manner an entire grid (centered about the input or computed firing 
attitude) can be considered for each ignition time. 

6. The user can elect to have the satellite post maneuver state vectors 
for each combination of ignition time and firing attitude (defined in 
items 1 through 5 above) stored in a sequential file. This file can 
then be accessed by other programs to perform various studies or 
to serve as a starting point for orbit determination after the apogee 
motor burn. 

The IMPMOT Program can use any one of four sources of satellite ephemeris 
information. 

1. FDS Orbit File 

2. ORBl tape 

3. EPHEM tape 

4. Card input for use in the Brouwer propagator (BPACKC) 

Any one of three sources of attitude information (if the user elects to input 
the firing attitude) can be used; 

1. FDS Attitude File 

2. SPINAT attitude prediction tape 

3. Card input 
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Either of two sources of attitude uncertainty information can be used 


1. FDS Attitude File 

2. Card input 
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SECTION 2 


INPUT DEFINITIONS AND FORMATS 

This section discusses the input required on FORTRAN logical unit number 5. 
Seventeen data cards are required to exercise fully all of the options of this 
program. The content and format of these cards are described in Tables 1 
through 12. Some additional comments are provided below. 

Cards lA, 2A, 3A, 4A, 5A, 7A and 8A are each data option indicator cards 
and each must precede its corresponding data card IB, 2B, 3B, 4B, 5B, 7B 
and 8B, respectively. If a particular B-card is omitted, then the correspond- 
ing A-card must also be omitted. 

Cards 2A and 2B are not applicable for the current version of the IMPMOT 
Program. These cards are reserved for future use and will be ignored by the 
program if they appear in the data stream. 

There are two circumstances when it is desirable to omit one or more data 
cards. 

1. If the scan upon apogee motor ignition time is desired, both cards 7A 
and 7B must be omitted. 

2. If the program is to compute the firing attitude for each ignition time 
(based upon either the circular orbit guidance scheme or the initial 
velocity scheme), both cards 8A and 8B must be omitted. 
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Card Group 1; FDS Interfaces 


2. 1 

Card 1 contains the FDS interface information. The parameters on this card 
(see Table 1) control the levels of the orbit and attitude information used, as 
well as the level of the FDS Manager Summary Report gehsrated. The levels 
referred to herein relate to the various entry points of a direct access 
sequential data set. 

2. 1. 1 Orbit Data Source Indicators 

For IORB2 = 0, either epoch elements (card input, an ORB-1 tape, or an 
EPHEM tape must be provided. 

For IORB2 = 1, satellite ephemeris information will be retrieved from the 
FDS Orbit File residing on FORTRAN logical unit number 20. Level LVLOBT 
of this file will be used. 

2. 1. 2 Attitude Data Source Indicators 

For IATT2 = 0, attitude information (if desired) is supplied either via Card 
8B or the SPINAT tape. 

For IATT2 = 1, satellite attitude information will be retrieved from the FDS 
Attitude File residing on FORTRAN logical unit number 25. Level LVLATT 
of this data set will be used. Whenever the firing attitude is input (either via 
card 8B or the Attitude File) , this input firing attitude will be used for all 
apogee motor ignition times. 
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2. 1. 3 FDS 'Manager Summary Report Indicators 


For LVLFDR = 0, the FDS Manager Summary Report will not be generated. 
For LVLFDR > 0, level, LVLFDR, of the FDS Manager Summary Report 
(see Table i) will be written to FORTRAN logical unit number 50. This data 
set is sequential direct access. No more than 100 lines, containing 72 alpha- 
meric characters each, will be written per level. A scratch sequential data 
set, on FORTRAN logical unit number 10, must be provided if LVLFDR ^ 0. 
For LVLFDR = -1, the next available level of the FDS Manager Summary 
Report will be written. For this option, the first eight bytes of each logical 
record on the data set defined by FORTRAN logical unit number 50 must be 
preihitialized with "NO DATA ". 

The parameter ISEXST controls the generation of the sequential excursion 
file. This file contains the satellite state vectors after the apogee motor 
bum for each combination of ignition time and firing attitude examined. For 
ISEXST = 0, the sequential excursion file will not be generated. 

For ISEXST = n, the sequential excursion file will be written on FORTRAN 
logical unit number 11. The program will assign reference numbers to each 
ignition time-firing attitude combination starting with "n". Each logical 
record (which corresponds to one ignition time-firing attitude combination) of 
this file contains the following information. 

1. Ignition time 

2. Burn-out time 
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3. Firing attitude 


4. Satellite position and velocity vectors after the apogee motor burn 

5. The osculating orbital elements after the apogee motor bum 

When specifying ISEXST = n, at least n-1 logical records must have been pre- 
viously written on the data set defined by FORTRAN logical unit number 11. 

The term "File 11" is synonymous with the sequential excursion file. 

2. 2 Card Groups lA and IB; Thrust Profile 

Card lA contains "99991" in card columns 1-5 (see Table 2). Upon reading 
this card, the program branches to the logic necessary to read Card IB. 

Card IB contains the apogee motor bum integration control parameters as 
well as the polynomial coefficients for the multisegment thrust profile. Card 
IB is physically composed of several data cards. The number of cards is 
determined by the input parameters which control the number of thrust seg- 
ments and the number of coefficients per thrust segment (parameters NSEG 
and NCOEFF, respectively). 

When reading card IB, the program begins by reading the number of thrust 
segments (NSEG) and the number of coefficients per each segment (NCOEF(i) 
where i = 1, 2, ... NSEG). A maximum of 10 segments is permitted. Each 
segment may have a maximum of 10 coefficients describing it. The polynomial, 
, which describes the thrust in the i**’ segment is defined by: 



where 


n — is the number of coefficients describing the i*'’ segment (n = 
NCOEF(i)). 

tj — is the segment time (details of this segment time will be given in 
subsequent paragraphs). 

From the above equation it is seen that the order of the polynomial describing 
the segment is NCOEF(i)-l. 

The program then reads a card containing the segment number (ISEG), the 
begin time and end time for the segment (BEGTM and ENDTM, respectively) 
and the integration time step to be used for the segment (TENCS).* Next, the 
coefficients describing the ISEG segment are read. The program will read in 
fields of 10 until NCOEF (ISEG) coefficients have been entered. 

The program will then read parameters ISEG, BEGTM, ENDTM, TINGS and 
COEFF (in the order described in the preceding paragraph) until the data for 
all NSEG segments have been entered. 

The segment time, t^, in equation 2. 2-1 is measured from the begin time for 
each segment. Thus, tj at the begin time of each segment is always zero. 

The relationship between the segment time for the i**' segment and the time 


*The begin times and end times for each segment are referenced to the 
beginning time of the apogee motor bum. Hence, the begin time for the first 

segment is 0. The beginning time of the i*** segment must equal to the end 

time of the i - 1®‘ segment. 
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measured from the beginning of the apogee motor burn, t, is given below 

t = tj + BEGTM (i) 

The program internally accounts for this time shift. The user must only 
remember that: 

1. The begin and end times are measured relative to the beginning of the 
apogee motor burn 

2. The coefficients describing the thrust polynomials (Cj j) are based 
upon segment time (tj). 

2. 3 Card Groups 2A and 2B: Reserve 

Cards 2A and 2B are not used in the current version of IMPMOT. Their 
presence is reserved for future improvements to this program. These cards 
will be ignored if input. 

2. 4 Card Groups 3A and 3B: Attitude and Guidance Data Card 

Card 3A contains "99993" in card columns 1-5 (see Table 2). Upon reading 
this card, the program branches to the logic necessary to read Card 3B (see 
Table 5). Card 3B contains spin axis uncertainty parameters (DELA and 
DEED); ignition time scan increment (DELTIM); attitude excursion parameters 
(SIGALP, SIGDEL, NSIGAL and NSIGDE); guidance option switch (LET); and 
indicators for attitude and attitude uncertainty sources (JATT and JERR).. As 
mentioned in Section 1, the user has one of two sources available for entering 
attitude uncertainty (either via the EDS Attitude File or via card input). 
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If JERR = 0, attitude uncertainty (in terms of the uncertainty in spin axis 
right ascension and declination) is entered via parameters DELA and DEED. 

If JERR 1 and the Attitude File option is selected (IATT2 = 1 on Card 1), 
the attitude uncertainty will be retrieved from level LVLATT (see Card 1) of . 
the FDS Attitude File. 

Through the parameters NSIGAL and NSIGDE, the user selects the range of 
the attitude excursion to be performed for each ignition time. If NSIGAL and 
NSIGDE are both zero, then the attitude excursion will not be performed. 

For NSIGAL 0 and NSIGDE ^ 0 the attitude excursion will be performed for 
each combination of: 

- NSIGAL « SIGALP<a<ap + NSIGAL • SIGALP 
8 j 3 - NSIGDE • SIGDEL < 5 < 5p + NSIGDE • SIGDEL 

where 

—is the input or computed spin axis right ascension 
—is the input or computed spin axis declination 
SIGALP and SIGDEL are defined in Table 5. 

For NSIGAL 4 0 and NSIGDE = 0, the excursion is performed for: 

- NSIGAL • SIGALP <a + NSIGAL • SIGALP 

5=6d 

For NSIGAL = 0 and NSIGDE 0, the excursion is performed for: 

a = «D 

5p - NSIGDE • SIGDEL < 6 < 6^ + NSIGDE • SIGDEL 
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The increments for the attitude excursions described above are: 

Aa = SIGALP 
A5 = SIGDEL 

If the user elects to input the firing attitude via the FDS Attitude File, the 
parameter JATT is set equal to 1 (whenever JATT = 1, the Attitude File 
option, IATT2, must be set equal to 1, see Card 1). 

The parameter LET indicates the guidance scheme to be employed. For 
LET = 0, the circular orbit guidance scheme will be used. The program will 
attempt to achieve a circular orbit for each apogee motor ignition time by 
trading off orbital inclination and node for eccentricity. The scan upon 
apogee motor ignition time will be terminated whenever the delta-V supplied 
by the apogee motor is insufficient to achieve a circular orbit (hence, the 
parameter DELSRM is ignored). Greater detail concerning this scheme can 
be found in the RAEMOT Program documentation in Reference 1. 

For LET = 1, the initial velocity guidance technique is employed. For each 
ignition time, the program will apply the apogee motor delta-V in the direction 
defined by the instantaneous velocity vector at the beginning of the burn. The 
range of the scan upon ignition time is controlled by the input parameter 
DELSRM (See Table 5). Greater detail concerning this guidance scheme can 
be found in Appendix A. 
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Card Groups 4A and 4B; Apogee Motor Characteristics 


2. 5 

Card 4A contains "99994" in card columns 1-5. Upon reading this card, the 
program will branch to the logic necessary for reading Card 4B. Card 4B 
contains the target orbit pararrieters and the apogee boost motor characteris- 
tics. The target orbit semimajor axis, eccentricity and inclination are used 
only for the purpose of computing the- residuals between the target orbit and 
the orbit resulting from the apogee boost motor burn and do not in any way 
affect the guidance scheme. The apogee boost motor impulse, XIMPLS, is 
defined by: 

XIMPLS = WEXP • I 


where: 

WEXP is the fuel weight expended during the apogee boost motor burn. I is 

oir 

the apogee boost motor specific impulse. 

An alternate expression for XIMPLS if thrust is assumed to be variable during, 
the apogee motor burn is given below 


XIMPLS = 



In (Wj/wj.) 


where 

T = thrust 

tg = apogee motor bum time 
w = weight flowrate 
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w = instantaneous spacecraft weight 


Wj = spacecraft weight at the beginning of the apogee motor burn 
Wj. = spacecraft weight at the end of the apogee motor burn 
2. 6 Card Groups 5A and 5B: Ground Station Data 

Cards 5A and 5B are optional. The default ground station data presented in 
Table 13 are assumed if these cards are not input. 

It is necessary that at least one ground station see the satellite at the apogee 
motor ignition time. If none of the ground stations can see the satellite, the 
ignition time is unacceptable and the program will proceed to the next card- 
stacked case (see Section 2. 13). 

Card 5A contains "99995" in card columns 1-5. Upon reading card 5A, the 
program will branch to the logic necessary to read Card 5B. Card 5B contains 
the ground station information (see Table 7). Up to 15 type 5B cards may 
follow a 5A card. Card 5A is read only at the beginning of this series. After 
the last 5B card in the series, a blank card is inserted. 

Station longitude is measured east from Greenwich from 0 to 360 degrees. 
Hence, a station at longitude 90 degrees west must be input as 270 degrees. 

If the station lies south of the equator, all non-zero latitude variables (LATD, 
LATM and LATS) must be input negative. 
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2. 7 Card Group 6A; Program Pause 

Card 6A is optional. It contains "99996" in card columns 1-5. Upon reading 
this card, the program enters a "pause" state and displays the number "7117" 
to the computer operator. 

To terminate the "pause" state, the computer operator enters any one-digit- 
character (non blank) via the computer console. During the "pause" state, 
the allocation of various I/O devices can be modified (e. g. different ORBl 
and EPHEM tapes can be mounted or dismounted within a job step). This 
option should only be used by an experienced user running the program in a 
"hands-on" mode. 

2. 8 Card Groups 7A and 7B; Firing Time 

Cards 7A and 7B are optional input cards. They are input only if the user 
does not want the program to scan upon apogee boost motor ignition time. 

Card 7A contains "99997" in card columns 1-5 (see Table 2). Upon reading 
Card 7A, the program branches to the logic necessary to read Card 7B. 

Card 7B contains the calendar date and time of the desired apogee motor firing 
time (see Table 8). 

The program assumes that the firing time input on this card is the desired 
burn-out time, tg^, for the apogee motor. The apogee motor ignition time 
(tf) corresponding to tgp is defined by: 

tf = tgo - TBOOST - BRNDLY 
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where 


TBOOST — is the apogee motor burn time (see Table 6) 

BRNDLY - is the apogee motor bum delay (see Table 6) 

2. 9 Card Groups 8A and 8B: Firing Attitude 

Cards 8A and 8B are optional input cards. They are input only if the user 
wishes to suppress internal calculation of apogee motor firing attitude and 
specify a particular firing attitude. Card 8A contains ”99998" in card columns 
1-5. 

Upon reading Card 8A, the program branches to the logic necessary to read 
Card 8B. Card 8B contains input attitude information (see Table 9). If the 
Attitude File option is used (lATT = 1, on Card 1 and the parameter JATT = 1 
(see Table 5), the variables ALPHA and DELTA on this card may be left 
blank. 

2. 10 Card Group 9; Apogee Search Interval 

Unless the desired apogee motor firing time is input (see Section 2. 8), the 
program must perform a search in order to find apogee of the transfer tra- 
jectory. The details of this search procedure are discussed in Section A. 1 of 
Appendix A. 

Card 9 contains the start and end times for the apogee search (see Table 10). 
The program begins its search for apogee of the transfer trajectory at the 
time defined by IBGRUN. If apogee does not lie within the time interval 
bounded by IBGRUN and lEDRUN, the search is terminated and the program 


19 



proceeds to the next card-stacked case (see Section 2. 13). By proper selec- 
tion of the start and end times, various apogees of the transfer trajectory 
may be examined. 

If the FDS Orbit File, ORB-1 tape, or EPHEM tape is used, care must be 
taken to ensure that neither IBGRUN or EEDRUN extend beyond the bounds of 
the available satellite ephemeris data. 

If the desired firing time is input, the input begin and end times will be 
ignored by the program. 

2. 11 Card Group 10: Transfer Orbit Characteristics 

Card 10 (see Table 11) contains the transfer orbit characteristics as well as 

the following indicators specifying: 

1. Transfer orbit elements tape 

2. Use of a spin axis prediction tape 

3. Perturbations to be considered for the Brouwer propagator 

4. Use of debug print 

5. Time reduction factor to be used during the apogee search 

If either the FDS Orbit File, ORB-1 tape or EPHEM tape is used, the follow- 
ing parameters may be left blank: TRSMA, TRECC, TRINC, TRMEAN, 
TROMEG, TRNODE and IPSET. 

2. 11. 1 Orbit Information Sources 

Card 10 in conjunction with Card 1 informs the program of the source of the 
satellite ephemeris information to be used. As mentioned in Section 1, the 
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IMPMOT Program has a choice of one of four sources for this informa- 
tion, 

1. FDS Orbit File 

2. ORBl tape 

3. EPHEM tape 

4. Card input for use by the Brouwer propagator (subroutine BPACKC) 
The option switches which flag each of these sources are described below. 

2. 11. 1. 1 Flags for the FDS Orbit File 

To use the FDS Orbit File, the user inputs IORB2 = 1 on Card 1 and indicates 
the level of the Orbit File to be used via the parameter LVLOBT (see Table 1). 
In addition the parameter TRSMA (input on Card 10) must be left blank. 

2. 11. 1. 2 Flags for an ORBl Tape 

An ORBl tape may be supplied in either a seven-track or nine-track format. 
Regardless of the format, the parameters IORB2 must be zero; TRSMA must 
be left blank; and lEPHEM is set equal to 0. 

If a seven-track ORBl tape is used, ITRACK is set equal to 7 (when using any 
other option, including either the Orbit File or EPHEM tape, ITRACK should 
be set equal to 9). The seven-track double-precision ORBl tape is supplied 
on DD name GENTAP. The seven-track tape will be converted to nine-track 
single-precision format and stored on FORTRAN logical unit number 17. If 
the debug print control (IBUG) is set, the conversion of the first 50 records of 
the ORBl tape will be printed. 
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If a nine-track ORBl tape is used, it must be in single-precision format. The 
nine-track single-precision ORBl tape is entered directly on FORTRAN logical 
unit number 17. The parameter ITRACK must also be set equal to 9. 

2. 11. 1. 3 Flags for an EPHEM Tape 

If an EPHEM tape is used, the parameters IORB2 must be zero; TRSMA must 
be left blank; and lEPHEM is set equal to 1. The nine-track single-precision 
EPHEM tape is then entered on FORTRAN logical unit number 17. 

2. 11. 1. 4 Flags for the Brouwer Propagator 

To input orbital elements via cards and use the Brouwer propagator, IORB2 
on Card 1 must be set equal to 0. The orbital elements are then entered via 
Card 10 via parameters TRSMA, TRECC, TRINC, TRMEAN, TROMEG and 
TRNODE. These input elements may b© either: 

1. Brouwer mean elements (IE LEM = 0 for this case) 

2. Osculating elements (IE LEM = 1 for this case) 

2. 11. 2 Attitude Information Source 

Card 10 in conjunction with Cards 1 and 3B inform the program of the source 
of input firing attitude data (see Sections 1, 2. 4 and 2. 9 for the implications 
of inputting the firing attitude). As mentioned in Section 1, the IMPMOT 
Program can use any one of three sources for this information. 

1. FDS Attitude File 

2. SPINAT attitude prediction tape 

3. Card input 
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The option switches which flag each of these sources are described 
below. 

2. 11. 2. 1 Flags for the FDS Attitude File 

In order to use the FDS Attitude File, the user inputs IATT2 = 1 on Card 1 
and indicates the level of the Attitude File to be used via the parameter LVLATT 
(see Table 1). The parameter ISPIN on Card 10 (see Table 11) must be set 
to zero, and JATT on Card 3B is set to one. In addition. Cards 8A and 8B 
are omitted. 

2. 11. 2. 2 Flags for the SPENAT Tape 

To use the SPINAT attitude prediction tape, the user inputs: IATT2 = 0 on 
Card 1; JATT = 0 on Card 3B; and ISPIN = 1 on Card 10. Cards 8A and 8B 
are omitted. The SPINAT tape is entered via FORTRAN logical unit number 
36. 

2.11.2.3 Flags for Card Input Attitude Data 

The user can elect to input the firing attitude via Cards 8A and 8B. Section 
2. 9 describes the format of these cards. To use this option the user inputs: 
JATT2 = 0 on Card 1; JATT = 0 on Card 3B; and ISPIN = 0 on Card 10. 

2. 12 Card Group 11: Transfer Trajectory Epoch 
Card 11 contains the epoch time of the transfer trajectory as well as user 
supplied comments (see Table 12). If the Orbit File, ORBl or EPHEM option 
is used, the epoch times may be left blank. 
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2. 13 


Multiple Card-Stacked Cases 


This program can, run any number of card-stacked cases per job step. The 
number of cases need not be known in advance. In order to run multiple cases, 
cards lA through 11 are repeated in sets for each case. 

If the Orbit File or Attitude File option is selected, their levels may not be 
altered among the stacked cases. After the desired number of cases have 
been stacked (which may consist of only a single case) , a card containing 
"99999" in card columns 1-5 is inserted. Figure 1 illustrates the data card 
stacking sequence for multiple cases. 
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Table 1. Card Group 1: FDS Interfaces (Sheet 1 of 2) 
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If LVLFDR 0, a scratch sequential data set on FORTRAN logical unit number 10 must be provided. 
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Table 2. Card Groups lA through 8A: Data Option Indicators 
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CARD 1 B is physically composed of several cards (note the indicated card columns). 




Table 2. Card Group IB: Thrust Profile (Sheet 2 of 2) 
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Table 5. Card Group 3B: Attitude and Guidance Data (Sheet 1 of 3) 
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If the Attitude File option is selected (IATT2 = 1), and JERR = 1, these parameters may be left blank. 
This parameter is recognized only if LET = 1. 



Table 5. Card Group 3B: Attitude and Guidance Data (Sheet 2 of 3) 
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If the Attitude File option is not selected, these parameters must be left blank. 




Table 5. Card Group 3B: Attitude and Guidance Data (Sheet 3 of 3) 



33 





34 




Table 7. Card Group 5B: Ground Station Data 
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Table 8. Card Group 7B: Firing Time 
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Table 9. Card Group 8B: Firing Attitude 
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Table 11. Card Group 10: Transfer Trajectory Characteristics (Sheet 1 of 3) 
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If either the Orbit File, ORB1 or EPHEM option is used, these parameters are left blank. 




Table 11 . Card Group 10 : Transfer Trajectory Characteristics (Sheet 2 of 3 ) 
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*lf either the Orbit File, 0RB1 or EPHEM option is used, these parameters are left blank. 

**The terms K2 and K4 are Brouwer eouivalents of the earth perturbation terms Jj and (see Reference 2 ). 
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Table 13. Default Ground Station Table 


STATION 

NAME 

NORTH 

LATITUDE 

EAST 

LONGITUDE 

ALTITUDE* 

M 

MINIMUM 

ELEVATION 

DEG 

MIN 

SEC 

DEG 

MIN 

SEC 

DEG 

ROSMAN 

35 

11 

46 

277 

7 

27 

876.00 

10.00 

FTMYRS 

26 

32 

54 

278 

8 

5 

9.00 

10.00 

MOJAVE 

35 

14 

49 

243 

6 

1 

921.00 

10.00 

ALASKA 

64 

52 

19 

212 

9 

40 

189.00 

10.00 

ORORAL 

-35 

-37 

-38 

148 

57 

11 

947.00 

10.00 

MADGAR 

-19 

0 

-25 

47 

18 

0 

1361.00 

10.00 

JOBURG 

-25 

-53 

-1 

27 

42 

28 

1 565.00 

10.00 

WNKFLD 

51 

26 

45 

359 

18 

14 

87.00 

10.00 

SNTAGO 

-33 

-8 

-56 

289 

19 

53 

681.00 

10.00 

LIMAPU 

-11 

-46 

-35 

282- 

50 

59 

34.00 

10.00 

QUITOE 

0 

-37 

-21 

281 

25 

16 

3578.00 

10.00 

NEWFLD 

47 

44 

29 

307 

16 

47 

112.00 

10.00 


'Measured above the mean equitorlal radius. 
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Figure 1 . Data-Card Stacking Sequence 
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SECTION 3 


OUTPUT DEFINITIONS 

This section discusses the printed output which has been subdivided into the 
following seven categories. 

1. Initial conditions and apogee motor characteristics 

2. Debug print (this will appear only if the debug print option is 
selected, see Card 10) 

3. Ignition-time "case" results* 

4. "Special case" summary 

5. FDS Manager Summary Report (this will appear only if ICOPY = 1 
and LVLFDR 7^ 0, see Card 1) 

6. Ignition time scan and firing attitude excursion table (this will appear 
only if an attitude excursion is performed, see Card 3B) 

7. Ignition time scan and attitude uncertainty table 

The printed output from categories 1 through 5 are directed to FORTRAN 
logical unit number 6. The output from categories 6 and 7 are directed to 
FORTRAN logical unit numbers 8 and 9, respectively. Table 14 presents an 
alphabetical list of each symbol, along with its definition and physical units, 


*The ignition-time "case" mentioned above is not to be confused with the card- 
stacked case mentioned in Section 2. 13. Within each card-stacked case, a 
scan upon apogee motor ignition time can be performed (see Section 1). Each 
ignition time which is examined within a particular card- stacked case is 
defined to be an ignition-time "case". 
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which appear in output categories 1, 2, 3, 4, 6 and 7. A "1" to the left of the 


symbol name indicates that it appears only if the debug print option is 
selected. Table 15 presents an alphabetical list of symbols appearing in the 
FDS Manager Summary Report (output category 5). 

3. 1 Initial Conditions and Apogee Motor Characteristics 

Figure 2 presents a sample of this portion of the output. Contained in this 

section are: 

1. Input begin and end times as well as user supplied comments (see 
Sections 2. 10 and 2. 12) 

2. Transfer orbit characteristics 

3. Target orbit characteristics 

4. Apogee motor characteristics 

5. Attitude uncertainties 

6. Thrust profile parameters 

7. Ground station data 

8. Source of the satellite ephemeris information in the transfer trajec- 
tory 

3. 2 Debug Printout 

The debug print precedes the ignition-time "case" with which it is associated. 
A sample of this debug print is presented in Figure 3. The parameters used 
in this portion of the output are noted by a "1" on Table 14. The debug print 
contains: 
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1. Inertial to orbital frame transformation matrices* 

2. Satellite position and velocity vectors prior to the apogee motor burn 

3. Delta-V vector supplied by the apogee motor** 

4. Apogee motor firing attitude 

5. Sensitivity of resulting orbit eccentricity to attitude errors and the 
residual semimajor axis and inclination (target minus achieved) 

6. Orbit true anomaly before and after the apogee motor burn 

If the firing attitude is input, either via card input, the SPINAT tape, or the 
FDS Attitude File, the variables DVB and DVBIJK (see Figure 3 and Table 14) 
will not be the delta-V vector applied in the inertial and orbital frames. This 
occurs because these variables are computed and output prior to encountering 
logic which imposes the fixed input attitude constraint. Whenever the firing 
attitude is input, these parameters should be ignored. 

*The orbital frame mentioned above is the local-vertical system defined by the 
satellite position and velocity vectors. The axes of this local-vertical system 
(f) '^and'i^) are defined by: 

-> 

(local vertical) 

IRI 

-*■ — > 

R X V 

^ (orbit normal) 

IR X VI 

f =^x1c (local horizontal) 

— y — ► 

where R and V are satellite position and velocity vectors, respectively. 

**The delta-V vector printed is based upon the magnitude of the delta-V as 
computed from the parameters input on Card 4B (see Section 2. 5 and Table 6). 



47 



3.3 


Ignition- Time "Case" Results 


Each card-stacked case (as mentioned in Section 2. 13) generates one or more 
ignition-time "cases". Each apogee motor ignition time which is examined 
within a card-stacked case constitutes an ignition-time "case". A sample of 
an ignition-time "case" output is provided in Figure 4. The information pre- 
sented for each "case" result includes: 

1. Ignition time case number and apogee motor ignition time 

2. Spin axis attitude at ignition 

3. Satellite position and velocity vectors prior to the apogee motor bum 
as well as the delta-V required to achieve a circular orbit at the 
ignition time 

4. Delta-V vector applied by the apogee boost motor 

5. The angle between the position and velocity vectors at the ignition; 
time 

6. Resulting orbit parameters 

7. The residual semimajor axis and inclination (target minus achieved) 

8. Sensitivity of resulting orbit eccentricity to attitude errors 

9. Ground station elevation angles at the apogee motor ignition time 

3. 4 Special "Case" Summary 

A "special case" summary is generated at the end of each card-stacked case. 
The case numbers which appear in this summary (see Figure 5) refer to the 
ignition-time "cases". The data in this output portion is self-explanatory. 
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FDS Manager Summary Report 


3. 5 

If ICOPY = 1 (see Card 1), a printed copy of the FDS Manager Summary 
Report is provided. A sample of this report is presented in Figure 6. Table 
15 presents the symbols and definitions used. Included in the FDS Manager 
Summary Report are: 

1. Orbit level used (OLVL) (if the Orbit File option is not used, this 
indicator will be zero) 

2. Attitude level used (ALVL) (if the Attitude File option is not used, 
this indicator will be zero) 

3. The orbit level corresponding to the attitude level used (AOLV) (if 
the Attitude File option is not used, this indicator will be zero) 

4. Total number of cases processed 

5. Epoch orbit elements 

6. Input spin axis attitude and uncertainties* 

7. Ignition time (day, hour, minute, and second in packed integer for- 
mat) 

8. Spin axis right ascension and declination 

9. Resulting orbit: 

a. semimajor axis 

b. inclination 

*If the spin axis is not input, the variables AW and DW are set equal to 180 
and 0, respectively. 
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c. right ascension of the ascending node 

d. perigee radius 

e. apogee radius 

3. 6 Ignition Time Scan Tables 

For each apogee motor ignition time examined, two tables will be generated; 

1. Attitude excursion table 

2. Attitude uncertainty table 
Each of these tables contains: 

1. Apogee boost motor ignition time 

2. Spin axis right ascension and declination at ignition time 

3. Resulting orbit: 

a. semimajor axis 

b. eccentricity 

c. inclination 

d. mean anomaly 

e. argument of perigee 

f. right ascension of the ascending node 

g. radius of perigee 

h. radius of apogee 

4. Reference number 

The term "FILE 11” which precedes each of these tables is synonymous with 
the sequential excursion file. 
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3. 6. 1 


Attitude Excursion Table 


If either NSIGAL or NSIGDE or both NSIGAL and NSIGDE are nonzero, an 
excursion upon firing attitude will be performed (see sample output in Figure 
7). The limits and increments for this excursion have been presented in 
Section 2. 4. The S 5 rmbol "E" stands for the increment of either spin axis 
right ascension or declination used for the attitude excursion. 

3. 6. 2 Attitude Uncertainty Table 

Referring to Figure 8 it is seen that for each ignition time a table is generated 
for the following attitude combinations: 


1. 

«D 

and 

2. 

«D 

+ U and 5r, 

3. 

“d 

- and 

4. 

«D 

and 6p + Ug 

5. 

«D 

and - Ug 

where: 




is the input or calculated spin axis right ascension at ignition 
time 

is the input or calculated spin axis declination at ignition time 
is the uncertainty in spin axis right ascension obtained from either 
Card 3B (variable DELA) or the Attitude File 
Uj is the uncertainty in spin axis declination obtained from either 
Card 3B (variable DEED) or the Attitude File. 
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The symbol "EAW" stands for uncertainty in spin axis right ascension. 
The symbol "EDW" stands for uncertainty in spin axis declination. 
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Table 14. Output Variables (Sheet 1 of 3) 


SYMBOL 

DEFINITION 

UNITS 

A 

ORBIT SEMIMAJOR AXIS 

KILOMETERS 

ALPHA 

SPIN AXIS RIGHT ASCENSION 

DEGREES 

ALPHAW 

SPIN AXIS RIGHT ASCENSION 

DEGREES 

AW 

SPIN AXIS RIGHT ASCENSION 

DEGREES 

C 

CASE NUMBER 


DEL 

UNCERTAINTY IN SPIN AXIS RIGHT ASCENSION OR DECLINATION 

DEGREES 

DELTA 

SPIN AXIS DECLINATION 

DEGREES 

''delvr 

VELOCITY INCREMENT REQUIRED TO OBTAIN A CIRCULAR ORBIT 

KILOMETERS 
PER SECOND 

Idvb* 

COMPONENTS OF THE APOGEE BOOST MOTOR DELTA-V IN THE 
INERTIAL COORDINATE SYSTEM (X, Y, Z COMPONENTS IN 
ROW FORMAT) 

KILOMETERS 
PER SECOND 

'’dvbijk** 

COMPONENTS OF THE APOGEE BOOST MOTOR DELTA-V IN THE 
LOCAL VERTICAL COORDINATE SYSTEM (1, J, K COMPONENTS 
IN ROW FORMAT) 

KILOMETERS 
PER SECOND 

DW 

SPIN AXIS DECLINATION 

DEGREES 

E 

1 

ORBIT ECCENTRICITY 


EAW 

UNCERTAINTY IN SPIN AXIS RIGHT ASCENSION 

DEGREES 

ECC 

ORBIT ECCENTRICITY 


EDW 

UNCERTAINTY IN SPIN AXIS DECLINATION 

DEGREES 

1 

ORBIT INCLINATION 

DEGREES 

lAFMOD 

FIRING MODE INDICATOR 

= 1 IMPLIES DEFICIENT ENERGY MODE 
= 2 IMPLIES SUFFICIENT ENERGY MODE 
= 3 IMPLIES EXCESS ENERGY MODE 


lATTUD 

INPUT ATTITUDE INDICATOR 

= 0 FIRING ATTITUDE IS TO BE COMPUTED 
= 1 INPUT ATTITUDE IS TO BE USED 


IFIRE 

INPUT FIRING TIME INDICATOR 

= 0 FIRING TIME IS TO RECOMPUTED 
= 1 INPUT FIRING TIME IS TO BE USED 



^Appears only if debug print option is selected. 


*lf firing attitude is input, the values of this parameter should be ignored since they are calculated prior to the impo- 
sition of the fixed attitude constraint. 

• •The element I is the component of OE LTA-V normal to the position vector; element J is the component of DE LT A-V 

normal to the orbit plan; element K is the component of DELTA-V along the position vector. If firing attitude is 
input, the values of this parameter should be ignored. 
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Table 14. Output Variables (Sheet 2 of 3) 


SYMBOL 

DEFINITION 

UNITS 

bPIRST 

FIRING TIME SCAN INDICATOR. 

= 0 IMPLIES THAT THE INDICATED FIRING TIME IS AT 
THE APOGEE OF THE TRANSFER TRAJECTORY. 

= 1 IMPLIES THAT THE INDICATED FIRING TIME IS PRIOR 
TO APOGEE OF THE TRANSFER TRAJECTORY 
= 2 IMPLIES THAT THE INDICATED FIRING TIME IS AFTER 
APOGEE OF THE TRANSFER TRAJECTORY 


INCL 

ORBIT INCLINATION 

DEGREES 

IORB1 

ORB1 USAGE INDICATOR 

= 0 IMPLIES THAT AN ORB1 TAPE IS NOT USED 
= 1 IMPLIES THAT AN ORB1 TAPE IS USED 


bSERCH 

DUPLICATE FIRING TIME INDICATOR (SEE IFIRST) 


M 

ORBIT MEAN ANOMALY 

DEGREES 

M ANOM 

cIrBIT mean ANOMALY 

DEGREES 

N 

ORBIT NODE 

DEGREES 

1 NODE LA 

SENSITIVITY OF RESULTING ORBIT ECCENTRICITY TO NEGATIVE 
ERRORS IN SPIN AXIS RIGHT ASCENSION. 

(DEGREE)-^ 

’ngdeld 

SENSITIVITY OF RESULTING ORBIT ECCENTRICITY TO NEGATIVE 
ERRORS IN SPIN AXIS DECLINATION 

(DEGREE)-’ 

OMEGA 

ORBIT ARGUMENT OF PERIGEE 

DEGREES 

PERD 

TRANSFER TRAJECTORY PERIOD 

SECONDS 

RA 

RADIUS OF APOGEE 

KILOMETERS 

RP 

RADIUS OF PERIGEE 

KILOMETERS 

RANOOE 

RIGHT ASCENSION OF THE ASCENDING NODE 

DEGREES 

^RAORBF 

POSITION VECTOR, PRIOR TO THE APOGEE BOOST MOTOR BURN, 
IN THE LOCAL VERTICAL SYSTEM. 

KILOMETERS 

SMA 

SEMIMAJOR AXIS 

KILOMETERS 

’ST 

TIME FROM TRANSFER TRAJECTORY EPOCH 

SECONDS 

’tanomr 

ORBIT TRUE ANOMALY AFTER THE APOGEE BOOST MOTOR 
BURN 

DEGREES 

’tanomt 

ORBIT TRUE ANOMALY BEFORE THE APOGEE BOOST MOTOR 
BURN 

DEGREES 

Waorbf 

VELOCITY VECTOR, PRIOR TO THE APOGEE BOOST MOTOR 
BURN, IN THE LOCAL VE RTICAL SYSTEM. (THE FIRST ELEMENT 
REPRESENTS THE COMPONENT NORMAL TO THE POSITION 
VECTOR; THE SECOND ELEMENT REPRESENTS THE COMPONENT 
NORMAL TO THE ORBIT PLANE; THE THIRD ELEMENT REPRE- 
SENTS THE COMPONENT ALONG THE POSITION VECTOR.) 

KILOMETERS 
PER SECOND 


^Appears only if DEBUG print option is selected. 
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Table 14. Output Variables (Sheet 3 of 3) 



SYMBOL 

DEFINITION 

UNITS 


^VBOOST 

APOGEE BOOST MOTOR DELTA-V 

KILOMETERS 
PER SECOND 

- 

’vc 

CIRCULAR ORBIT VELOCITY 

KILOMETERS 
PER SECOND 

. 

VELAPO 

VELOCITY AT APOGEE OF THE TRANSFER TRAJECTORY 

KILOMETERS 
PER SECOND 


W 

ORBIT ARGUMENT OF PERIGEE 

DEGREES 


^Appears only if DEBUG print option is selected. 
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Table 15. FDS Manager Summary Report Variables 


SYMBOL 

DEFINITION 

UNITS 

A 

SPIN AXIS RIGHT ASCENSION AT IGNITION 

DEGREES 

ALVL 

ATTITUDE LEVEL (IF THE ATTITUDE FILE OPTIDN IS 
NOT SELECTED, THIS VARIABLE WILL BE ZERO) 


AOLV 

ORBIT LEVEL ASSOCIATED WITH ALVL (IF THE ATTITUDE FILE 
OPTION IS NOT SELECTED, THIS VARIABLE Wl LL BE ZERO). 


AW 

INPUT SPIN AXIS RIGHT ASCENSION* 

DEGREES 

D 

SPIN AXIS DECLINATION AT IGNITION 

DEGREES 

DW 

INPUT SPIN AXIS DECLINATION* 

DEGREES 

E 

ORBIT ECCENTRICITY 


EAW 

UNCERTAINTY IN SPIN AXIS RIGHT ASCENSION 

DEGREES 

EDW 

UNCERTAINTY IN SPIN AXIS DECLINATION 

DEGREES 

1 

ORBIT INCLINATION 

DEGREES 

M 

ORBIT MEAN ANOMALY 

DEGREES 

N 

ORBIT RIGHT ASCENSION OF THE ASCENDING NODE 

DEGREES 

OLVL 

ORBIT LEVEL (IF THE ORBIT FILE OPTION IS NOT SELECTED, 
THIS PARAMETER WILL BE ZERO) 


P 

PERIGEE RADIUS 

KILOMETERS 

RA 

APOGEE RADIUS 

KILOMETERS 

S 

ORBIT SEMIMAJOR AXIS 

KILOMETERS 

W 

ORBIT ARGUMENTOF PERIGEE 

DEGREES 


*lf the spin axis attitude is not input, AW and DW will appear as 180.000 and 0.0, respectively. 
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Figure 2. Sample Initial Conditions Printout 
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Figure 3. Sample Debug Printout for Ignition-Time "case 
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Figure 4. Sample Ignition-Time "Case" Results 
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Figure 5. Sample "Special Case" Summary 
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Figure 6. Sample FDS Manager Summary Report 



IGNITION TIME ATTITUDE SCAN ALPHAS* OELTAW A £ I M W N APOGEE PERIGEE REFNO 

9 2£ 12 £7 2 AM .DM 69.065 -22*^73 233914. 0.03064 29.230 179.902 142.433 20C.643 241060. 226746. 1 


< . i ■ . ! ! • I ■ ■ ■ ; . , 

(M<o<u)toscoooOtf>ONCo>o««Nn« 0 (>QM(Mn^to<oNOOi 4 «tf><osiDO>o«MCijONa 9 ‘o«i(\jn^(n 4 >o<fNn« 

p.z..i->. 4 ..«>N(UN(\i(MZ(\jn'nmnr)r>nnf 4 Z«« 4 « 44 ' 4 iAinu)Zinioio«<o«' 0 < 0 ' 0 <OZNNNK. 

tt. IL' tt.. *U.~ TtL*' 

(IJ ' tU III • • tli 


<OOtf)iO 4 >< 0 ( 0 < 0 <O**< 0 'O’ 0 tn< 0 <O< 0 U)i 0 O*"O' 4 * 0 IA' 0 ' 0 ' 0 tA'O 4 >*~<Ci 0 Oi 0 O 4 O(nO< 0 ***-< 0 < 0 O( 0 C< 0 < 0 tn< 0 < 0 *«O’C<Ol 0 

MC\J(MNMCMMC\J(Ma(M(MCMOJ(MMMM(VN(Z<VJ(>J<NjMNOiC\J<MM(VjQ:OJCgAIN<U(M(M<VJ(VjOi(X(MNCg(MrgN(VJA)(NJ(Vja:M(N|(VICg 

<\|<VJ<M(MNNCVJCM<MtlJOJN<UN(M<\J<M<Vi(VI(MUiM<\J<MMMM(MOJNNtUNN(UN<M(MMA)C\tNIUNAIN(\JN(MCMMCVI(\jUJ(V>(V)<MN 


for-foh'r)..w^N 

cutr>«^tf> 4 <ocMNin 

<m(Mr^<vtnoo><v 


<o<Mc^Ma)lurt M<« 
<sr)K<o(\)aiocjifO 
« 9 nr)r>« 4 )ou>o 


^r>h>v>N«cotOtf><o 

a)«atf>ov>aio<c<oo 


m HI o> H> o tf) m .. v* 10 

nw<uosiorar^p>« 

■iMAjoimtoinaiHiiv 


oair»<OMOv 9 iooi (viov*i^- 


<^ 0 . 4 V‘OimOO.« 

0 >OO< 7 tOOO<JO 

. 4 (MOI. 4 CMOiC\INN 


co.«(\l<Hr)Nao<j>o 

oin\ow'*>o<oo< 


OOOOOOOO 0-0 


<oin«r^roa 3 <\jO(Ma> 

n^oov«for> ^ *n •c 

%«U><M'O<UHfl 0 'llUO 


000 <j 000000 

N.h<M<M. 4 (MMN<MN 




000<^0(j0000 OOOOOOOOOQ oooo 

NtHMilDINNMCMN N.«(M(M(UCM<M(Uni(M NNMN 

90 ‘ 0 ' 000 j<t ^00 «h-(VIN 

oa>« 0 < 00 o^s> 0'0 inoa>Aj<otf)innp.N oiOnhi- 

<DOi*)aiO>(Dtf)( 0 <MN MO^(M(UMO«a) 0 > ui**Ai.< 


Oh.* 0 'u<\io>rMr-u 

fV},..«MO( 000 < 0 - 

«sOU>«*U'HIUCMCU 


CM (7> Ml 1^ <0 «H n «M HI M 

oinK«ifiuoo(r)r« 

<U \H « U V> (0 ri HI .4 


<H <0 .. 41 Ml (H « <0 HI 
Onim^OOlHli^CHOI 


oa>ni«*i<i>(oN ..«o 


«OvHir^(D«Oi<M>VH> 
OlfOiulO«(DNl|iUi(\] OOOO 

KMPiMii^r^uiaiMa) r^oi. 4 c\i 


Hiuiui 4 i<'‘r‘oO\^ M'V'iiuiair^r^'^O^ 


■ 1 ) u) HI r* M> VI o w> 


i<ciai>«r^i/>..r>*.Mi imhiui.. 
HI .1 HI HI P* r> O O (A P> HI ui ui 


Hj(M«noojtnfvjif) ..ovf«(vi.«p>cu<uin OHi<M^Mio«U)ior« 0 ‘i^ko.>oi9p>h)oi coohiiolpoimiionid 

(u.<^or)HioiHiiu O 40 '('i’ 0 o 4 <ooir» «o<viu>«.«^>r>H>ti> 0 (\i<o..Hiiui'>h‘..'0 Niouio<«p>..h>oio vo^nio 

vp«a>oN<OHiO« OP'M’Omo^fOaoftj .<a>r)r>if>>.«ininoi*) ocd<vjn<oou>«o>p> u>N«.u>u)v>«r)(00J (uavo 

(j>0'9)0‘a«ooo‘0 o<u(D^'OoNal^>^ oooa)o>«>a)i>a><D <s<O(0p>svNOr4>N oocpopoioocdo oodos 

OIMIMIVIHIMINPIM <V(M<MWOI(M0fCH<HM OlM<MCM(VMlVCM(M(M CVOJIMIHOICVICVJOIMOJ IHOlCMCM- 

0 - 400 «^ 0 > — lOO 4 N- 0 < 04040 '»* moi- — o«nN®o® «IOO*N«*-i®^S ♦ooto 


..SfflO^'OWT — n ^r».« 0 ‘ 04 <J >®'00 ^tMOOC^'OmOM <M 0 >«tf>-"N® ri® 0 > 0 ‘ 0 >'OiOiOK <0 lONOO* 

^ 00 'O^■•-NP'.-W .«<VJ. 4 SM*-P^OO’ 0 ' .-f* 0 <VS— NOO*" — (MOIOi.himcOO'®® — OI — C^J^.«^J^ 0 ©C^ -*.N<go 

CUOtfiHlUiliPUlViU* (PVI*OiUs;>V<HI<J'®® U> 90 >aiV>V>®®®HI UIUJUIV>® 9 > 0 HI®® 

ror’fOfO'onfonn nnnniofoionnri nromfonnn'ocoo nr>iop>fOior*ionro rinnionionroiojo lorofoio 

rofnnnroriM)nn<fnn"innM>iOPono<fnrtnnio»oioDior>^r>fOioionioiofOioio<iMr)P^'OP>(or>ro(oro<ioioior> 

MftJCMftiMNWNN CJN<M<NMtVI<MW<MN MMWWNNWMNM WNNfMMNf'iNNW NNCgCVl' 


in®inu)ioin'niniP oooooooooo cmmcmimcmcmnojcijn c^c^9'(^(7>Ol;•o^c^<^ (mn<mn 


LiUJUJtUUJUJUiUllUH uUJ(jmU>Uf(jlU(U‘ 

I I I U I I • < 


I UJ UJ III UJ ® z 

I VI .. .. O M < 

I I \J 


UJ UJ UJ UJ I 

.. .4 ^ O • 
< I I 


UlLJUJliJlillUUJUl ®3 (UUllUlUtiiUJW'HJU): 

..•H..OOO- 4 . 4 . 4 H' M. 4 . 4000 ......I 

III - I I I 


62 


Figure 7. Sample Ignition Time and Attitude Scan Table 
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Figure 8. Sample Ignition Time and Attitude Scan Table (attitude uncertainty only) 



SECTION 4 


JOB CONTROL LANGUAGE 

The Job Control Language (JCL) required to execute the IMPMOT Program on 
the IBM 360 series computer system at GSFC is presented in this section. 

For the purpose of discussion, it is assumed that an object file of the IMPMOT 
program resides on the first file of a nine-track, standard-labeled magnetic 
tape, YOURTAPE. 

The examples presented in the following subsections are the ones which the 
user is most likely to encounter. Additional combinations can be formed, 
using the examples as a guide, if the user keeps in mind the following rules 

1. DD names FTlOFOOl and FT50F001 are used for the FDS Manager 
Summary Report. FTlOFOOl defines a scratch sequential data set. 
FT50F001 identifies the FDS Manager Summary Report assigned to 
the IMPMOT Program. 

2. DD name FTllFOOl is reserved for the sequential excursion file. 

3. DD names FT17F001 and GENTAP are used in conjunction with 
ORBl and EPHEM tapes (see Section 2. 11. 1 for the internal flags 
which must be set). 

A. If a seven-track ORBl tape is used, it is entered via GENTAP. 
The program will convert the seven-track tape to nine-track 
format and store the result on FT17F001. 
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B. If a nine-track ORBl tape is used, it is entered directly via 
FT17F001. 

C. If a nine-track EPHEM tape is used, it is entered via FT17F001. 

4. DD name FT20F001 defines the FDS Orbit File. 

5. DD name FT25F001 defines the FDS Attitude File. 

6. DD name FT36F001 defines the SPINAT tape.* 

Seven examples are presented. 

1. JCL using the Orbit File and the Attitude File 

2. JCL using an ORBl tape and the Attitude File 

3. JCL using an EPHEM tape and the Attitude File 

4. JCL using the Brouwer propagator and the Attitude File 

5. JCL using the Orbit File and a SPINAT tape 

6. JCL using the Orbit File without the Attitude File 

7. JCL without the FDS Manager Summary Report 

In examples one through six, it is assumed that the FDS Manager Summary 
Report will be generated. For these examples it is further assumed that: 

1. The Attitude File data set, with data set name, ATTITUDE. FILE, 
resides on a disk. 

2. The Orbit File data set, with data set name, ORBIT. FILE, resides 
on a disk. 

*The SPINAT tape indicated above, refers to the attitude prediction tape men- 
tioned in Section 2. 11. 
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3. The SPINAT tape is a nine-track non-labeled magnetic tape with 
volume serial number SPINAT. 

4. The EPHEM tape is a nine-track non-labeled magnetic tape with 
volume serial number EPHEM. 

5. The EDS Manager Summary Report file, with data set name 

EDS. MANAGER. SUMMARY, is a preallocated data set residing on a 
disk. 

6. The Sequential Excursion Eile, with data set name, EXEILE, is a 
preallocated data set residing on a disk. This file is assumed to be 
generated for each example presented in this section. 

4. 1 JCL Using the Orbit and Attitude Eiles 

Card Ref. 
No. 

//LINK EXEC LINKGO, REGION. GO=300K 1 

//LINK.SYSLIN DD UNIT=2400-9 ,DISP= (OLD, KEEP) ,DSN=IMPMOT. OBJECT, 2 

// LABEL=(I,SL) ,DCB=(LRECL=80,RECFM=FB,BLKSIZE=3200) ,V0L=SER=Y0URTAPE 3 
//GO.FT05F001 DD DDNAME=DATA5 4 

//GO.FT06F001 DD SYSOUT=A,SPACE=(CYL, (3,2)) , 5 

// DCB=(LRECL=137, RECFM = VBA,BLKSrZE = 7265) 6 

//GO.FT08F001 DD SYSOUT=A,SPACE=(CYL, (3,2)) , 7 

// DCB=(LRECL=137,RECFM=VBA BLKSIZE=7265) 8 

//GO.FTlOFOOl DD UNIT=2314 ,DISP=(NEW, DELETE) ,DSN=&&SCRATCH, 9 

// SPACE=(CYL,(2,2)),DCB=(LRECL=72,RECFM=FB,BLKSIZE=7200) 10 

//GO.FTllFOOl DD UNIT=2314 ,DISP=SHR,DSN=EXFILE, 11 

// DCB=(LRLCL=164,RECFM=VBS,BLKSIZE=6564) ,VOL=SER=XXXXXX 12 

//GO.FT17F001 DD DUMMY 13 

//GO.FT20F001 DD UNIT=2314,DISP=SHR,DSN=ORBIT.FILE,VOL=SER=XXXXXX 14 

//GO.FT25F001 DD UNIT=2314 ,DISP=SHR,DSN=ATTITUDE .FILE,VOL=SER=XXXXXX 15 

//GO.FT36F001 DD DUMMY 16 

//GO.FT50F001 DD UNIT=2314,DISP=SHR,DSN=FDS. MANAGER. SUMMARY, 17 

// VOL=SER=XXXXXX 18 

//GO.DATA5 DD* 19 

//GO.GENTAP DD DUMMY 20 

/* ENTER DATA CARDS HERE 
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JCL Using an ORBl Tape and the Attitude File 


4. 2 

An ORBl tape may be in either a seven- or nine-track format. An example of 
each is presented below. Section 2. 11. 1. 2 describes the use of ORBl tapes. 

4. 2. 1 Seven-Track ORBl 

The changes to the JCL, as presented in Section 4. 1, needed to execute with a 
seven-track ORBl tape are as follows. 


Card Ref. 
No. 


//GO.FT17F001 DD UNIT=2314,DISP=(NEW, DELETE) ,DSN=&&NINE, 13 

// SPACE=(CYL, (5,2)) ,DCB=(LRECL=5608,RECFM=VS,BLKSIZE=5616) 13A 

// GO.FT20F001 DD DUMMY 14 

//GO.GENTAP DD UNIT=2400-7 ,DISP=(OLD,KEEP) ,DSN=SEVEN, 20 

// LAB£L=(1,BLP) ,DCB=(DEN=1,RECFM=V, BLKSIZE=4250) , VOL=SER=TAPE 20A 

//GO.FT20F001 


Where DD name GENTAP points to the seven-track double-precision ORBl 
tape. 

4. 2. 2 Nine-Track ORBl 

The changes to the JCL presented in Section 4. 1 which are required to execute 


with a nine-track single-precision ORBl tape are as follows 


Card Ref. 
No. 


//GO.FT17F001 DD UNIT=2400-9 ,DISP= (OLD, KEEP) ,DSN=NINE , 

/ / LABEL= ( , BLP) , DCB= ( LRECL=5 608 , RECFM=VS , BLKS IZE=5 616 ) , V0L= SER=TAPE 
//GO.FT20F001 DD DUMMY 


13 
13A 

14 


4. 3 JCL Using an EPHEM Tape and the Attitude File 

If an EPHEM tape is desired, the changes to the JCL presented in Section 4. 1 


are: 
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Card Ref. 
No. 


//GO.FT17F001 DD UNIT=2400-9,DISP=(OLD,KEEP) ,DSN=EPHEM, 13 

// LABEL=(1,BLP) ,DCB= (DEN=2 ,LRECL=2808 ,RECFM=VS , BLKSIZE=2816) , VOL=SER=EPHEM 13A 

//GO.FT20F001 DD DUMMY 14 


4. 4 JCL Using the Brouwer Propagator and the Attitude File 

The changes to the JCL presented in Section 4, 1 when using the Brouwer 


propagator are: 


Card Ref . 
No. 


//GO.FT20F001 DD DUMMY 14 


4. 5 JCL Using the Orbit File and SPINAT Tape 

The changes to the JCL presented in Section 4 . 1 when using the SPINAT tape 


in lieu of the Attitude File are as follows: 


Card Ref. 
No. 


//G0.FT25F001 DD DUMMY 15 

//GO.FT36F001 DD UNIT=2400-9 ,DISP= (OLD, KEEP), DSN=SPINAT, 16 

// LABEL=(1,BLP) , VOL=SER=SPINAT ,DCB= (LRECL=2808,RECFM=VS ,BLKSIZE=5616) 16A 


4. 6 JCL Using the Orbit File Without the Attitude File 

The change to the JCL presented in Section 4. 1 necessary to delete the Attitude 


File is presented below. 


Card Ref. 
No. 


//GO.FT25F001 DD DUMMY 


15 
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JCL Without the FDS Manager Summary Report 


4. 7 

The changes to the JCL in Section 4- 1 required if the FDS Manager Summary 
Report is not requested are: 


Card Ref. 
No. 


//GO.FTlOFOOl DD DUMMY 9 

//GO.FT50F001 DD DUMMY 17 
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APPENDIX A 


INITIAL VELOCITY GUIDANCE TECHNIQUE 

This appendix describes the logical sequence followed when the initial velocity 
guidance technique is employed. The major computational steps which the pro- 
gram follows are presented below. 

1. Compute the satellite state vector at apogee of the transfer trajectory 

2. Compute the satellite state vector at the ignition time 

3. Compute the apogee motor firing attitude 

4. Propagate the satellite state vector through the apogee motor burn* 

5. Perform the attitude excursion* 

6. Superimpose the attitude uncertainties upon the attitude computed in 
Step 3* 

7. Increment the ignition time 

8. Repeat Steps 2 through 7 until the ignition time scan is satisfied 
Whenever the desired firing time is input (see Section 2. 8), Step 1 is replaced 
by: 

lA. Compute the satellite state vector at the input firing time. In addition, 
Steps 7 and 8 are bypassed. 


*If the user requests that the sequential excursion file be generated, the 
satellite state vector (as defined by both the position and velocity vectors and 
the osculating orbital elements) are stored on this file at each of these steps. 
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A. 1 Satellite State Vector at Apogee 

The IMPMOT Program has four sources of satellite ephemeris information in 
the transfer trajectory 

1. FDS Orbit File 

2. ORBl tape 

3. EPHEM tape 

4. Card input for use in the Brouwer propagator 

Each of these sources provide the osculating orbital elements at epoch (a^, e^, 
nig, and 12^) as well as satellite position and velocity vectors as a function 
of time.** 

From and e^, the program makes an initial estimate of the satellite position 
and velocity at apogee (R^ and V^, respectively). 



The program begins its search for apogee of the transfer trajectory at the 
input begin time (see parameter IBGRUN on input Card 9). The satellite posi- 
tion and velocity vectors at the begin time are obtained and a comparison is 
made: 


**In the case of the FDS Orbit File and the Brouwer propagator, the satellite 
position and velocity vectors as a function of time involve additional calcula- 
tions but are nevertheless available. 


72 



!R[- I <200. kilometers 

|v| - Vg j < 0.2 kilometers per second 


where; 

R - is the satellite position vector at time t 
V - is the satellite velocity vector at time t. 

If both tests A-1 and A-2 are failed, the time, t, is incremented by: 

t = t + At 


where: 


A-1 


A-2 


A-3 


At - is the input time increment DELTIM (see Card 3B). 

The position and velocity vectors at this new time are obtained and tests A-1 
and A-2 are repeated. 

If either test A-1 or A-2 is passed, the satellite is in the vicinity of apogee 
and an additional test is made 


cos 


-1 


R-V 


IRI IVI 


> 7 , 


radians 


A-4 


If test A-4 is failed, equation A-3 is used to increment the time and the search 
is continued. If test A-4 is passed, the satellite has reached or passed 
apogee. If the equality portion of A-4 is met, apogee has been found and the 
search is terminated. If the inequality portion of A-4 is met, the satellite has 
passed apogee in which case one time increment is subtracted (i. e. t = t - At). 
The time increment is then reduced by the input parameter IDIVID (see Card 


10 ). 
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At = DELTIM/IDIVID 


A-5 


where 

IDIVID is the time reduction factor (see Card 10) 

The time is then incremented according to equation A-3 with the new time in- 
crement described in equation A-5. The search continues (i. e. R and V are 
obtained at the new time and comparisons A-1, A-2 and A-4 are made as in- 
dicated above) until A-4 is again satisfied. Once A-4 has been satisfied for 
the second time, with the reduced time increment, the apogee search is 
terminated. It is noted that the time increment was reduced only once during 
the apogee search scheme. 

Once the apogee search scheme has been satisfied, the satellite position and 
velocity vectors at apogee (R^ and V^, respectively) are available as well as 
the time at which the satellite reaches apogee of the transfer trajectory (t^ ). 
If the desired firing time is input, the apogee search logic described above is 
circumvented. The satellite position and velocity vectors at the desired igni- 
tion time (Rj. and Vj., respectively) are obtained at the requested firing time 
(tf). 

A. 2 Satellite State Vector at the Ignition Time 

The Program assumes that the apogee time, t^, is the desired burn-out time 
(if the firing time is input, the burn-out time is assumed to be tj.). The 
apogee motor ignition time, tjg, is computed from the burn-out time by: 

^ig “ ^Bo “ 


74 



where: 


tg is the apogee motor bum time (see parameter TBOOST on Card 
4B) 

tp is the burn delay time (see parameter BRNDLY on Card 4B) 
tg^ is the desired burn-out time (either or t^) 

The satellite position and velocity vectors (Rj and Vj, respectively) at time 
t.g are then obtained from one of the four sources indicated in Section A. 1. 

A. 3 Compute Firing Attitude 

The initial velocity guidance scheme applies the thrust vector in the direction 
defined by Vj. Hence a unit vector along the positive spin axis, is defined 
by (note that the thrust vector is assumed to be aligned along the negative spin 
axis) : 

A-7 




IV, 1 


where: 




The spin axis right ascension, a, and declination, 6, (which define the firing 
attitude) are computed by: 


= tan-ije^y/e^^| 

A-8 

6d = sin-1 

A-9 
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A. 4 Propagate the Satellite State Vector Through the Apogee Motor Burn 


The satellite state vector is propagated through the apogee motor burn by 
using an 8*'^ order Runga-Kutta technique to integrate the satellite equations of 
motion during the burn. (Reference 3 presents the satellite equations of motion) 
The satellite position and velocity vector at the beginning of the bum are R. and 
V^, respectively. The thrust vector, T, acting upon the satellite is defined by: 

T = -T-e^ A-10 

where: 

T - is generally a function of time and is computed based upon the 
input thrust coefficients (see Section 2. 2) 

A, 5 Perform the Attitude Excursion 

If the user requests the attitude excursion, each combination of spin axis right 
ascension, a, and declination, 5, in the following range will be taken for the 
apogee motor firing attitude. 

~ n • a < a < + n • a 

u Ql D oc 

5d ~ m • Oj < 6 < 5 d + • <^6 

where 

(Xp - is defined by equation A-8 
Sp- is defined by equation A-9 

cr - is the increment for the excursion upon right ascension (see 

a 

parameter SIGALP on Card 3B) 
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Og - is the increment for the excursion upon declination (see parameter 


SIGDEL on Card 3B) 


m and n are arbitrary integers (see parameters NSIGAL and NSIGDE 
on Card 3B) 

The unit vector along the spin axis is defined by: 

=> (cos a cos 5, sin a cos S, sin 6) A-11 

For each firing attitude in the range described above, the integration described 
in Section A. 4 is performed. 

A. 6 Examine Attitude Uncertainties 

The effects of attitude uncertainties upon the orbit resulting from the apogee 
motor burn are determined by superimposing these uncertainties upon the 
firing attitude computed by equations A-7, A-8 and A-9. Four cases are 
examined. 



, Of 


« = «D - 



a = «D 

. 5 = 5^ 


a = «D 

, 5 = 6^ 



where 


ttp - is defined by equation A-8 
5p -is defined by equation A-9 

-is the uncertainty in spin axis right ascension (see parameter 
DELA on Card 3B) 
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Ug - is the uncertainty in spin axis declination (see parameter DELD 
on Card 3B). 

For each firing attitude defined above, the integration described in Section A. 4 
is performed. 

A. 7 Increment Burn-Out Time 

The apogee motor bum-out time, is then incremented by an amount 
DELTIM (see Card 3B) and the sequence described in Sections A. 2 through 
A. 6 is repeated for each tg ^ in the range: 


^ ^BO 


<t, + At 


where 

t^ - is the time at which the satellite reaches apogee of the transfer 
trajectory (see Section A. 1) 

A tj-is the time period for the ignition time scan (see parameter 
DELSRM on Card 3B) 

A. 8 Logical Flowchart 

A descriptive flowchart illustrating the computational sequence described in 
Sections A. 1 through A. 7 is provided in Figures A-1 through, A-5. Each of 
these figures is outlined below. 

1. Figure A-1 illustrates the apogee search logic described in Section 

A. 1. 

» 2. Figure A-2 presents the ignition time and firing attitude computations 
presented in Sections A. 2, A. 3 and A. 4. 
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3. The attitude excursion logic described in Section A. 5 is illustrated in 
Figure A-3. 

4. Figure A-4 illustrates the attitude uncertainty logic described in 
Section A. 6. 

5. The apogee motor ignition time scan logic is illustrated in Figures 
A-4 and A-5. 
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Figure A-1 . Apogee Search 
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Figure A-2. Computed Ignition Time and Attitude 
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Figure A-3. Attitude Excursion (Sheet 2 of 2) 
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Figure A-4. AtHfude Uncerfainfy (Sheet 1 of 2) 
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Figure A-4. Attitude Uncertainty (Sheet 2 of 2) 
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